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Abstract 



By choosing l H nucleus in Carbon-13 labelled trichloroethylene as 
one qubit environment, and two 13 C nuclei as a two-qubit system, we 
have simulated quantum decoherence when the system lies in an entan- 
gled state using nuclear magnetic resonance (NMR). Decoupling tech- 
nique is used to trace over the environment degrees of freedom. Exper- 
imental results show agreements with the theoretical predictions. Our 
experiment scheme can be generalized to the case that environment is 
composed of multiple qubits. 
PACS number(s):03.67 

Quantum decoherence is a purely quantum-mechanical effect through 
which a system loses its coherence behavior by getting entangled with its 
environment degrees of freedom [1][2]. A real quantum system always in- 
teracts with its surrounding environment. The time evolution induced by 
the interaction introduces entanglement between the system and environ- 
ment when the system initially lies in a superposition of states. When the 
state of the system is described by a reduced density matrix through trac- 
ing over the environment degrees of freedom, decoherence makes the off- 
diagonal matrix elements approach 0, and leaves diagonal ones unaltered. 
Quantum decoherence is thought as a main obstacle for experimental imple- 
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mentations of quantum computation and has been widely studied in theory 
[3]- [5]. Some theoretical and experimental schemes have also been proposed 
to solve the problem of decoherence in order to fully use coherence in quan- 
tum information [6] -[8]. C.J. Myatt. et al observed and studied decoherence 
by coupling an ion in a Paul trap to a reservoir that can be controlled [9] [10]. 
D. G. Cory et al simulated decoherence by magnetic field gradients [7]. In 
this paper, we will simulate decoherence when the system lies in an entan- 
gled state. Although we choose one qubit to simulate the environment, our 
scheme can be generalized to the case that the environment is composed of 
multiple qubits. 

We choose carbon 13-labelled trichloroethylene (TCE) dissolved in d- 
chloroform as a sample. TCE's structure is shown in Fig.l. The two 13 C 
nuclei constitute the system, and the X H nucleus is viewed as environment. 
They represent qubits 1, 2, and 3, respectively. The influence of the other 
nuclei can be ignored. The Hamitonian of the system and environment is 
represented as 

H s+e = -coJl - uj 2 I 2 z - uj 3 I 3 z + 2tt J 12 l\ll + 2irJ 23 I 2 z I 3 z + 27rJ 13 / z 1 / 3 , (1) 

where Uk/2ir(k = 1, 2, 3) are the resonance frequencies of the 3 spins, I z are 
the matrices for the z-components of the angular momentum, and % is set 
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to 1. J 12 , for example, denotes the coupling constant between spins 1 and 
2. Before the system lies in an entangled state, 1 H nucleus is decoupled. J 23 
and J13 are averaged to zero. The system becomes a closed system whose 
Hamitonian is represented as 

H s = -uJl - uj 2 I 2 z + 2nJ l2 I l z I 2 z . (2) 

According to Ref. [11] and the experimental technique we used for heteronu- 
clear system [12], the pulse sequence [a] 2 .— [grad] z — [7r/4]^' 2 — 1/4 J i2 — M^' 2- 
1/4 J 12 — [— 5ir/6]y' 2 — [grad] z transforms the system from equilibrium state 
to the pseudo-pure state | |>i | |> 2 represented as p e f — 1\ + 1% — 2l\l 2 - 
Here [a] 2 , for example, refers to a radio- frequency (rf) pulse on spin 2, ori- 
ented along x-axis. The evolution caused by the pulse is denoted as e iaI *. 
[grad] z denotes a gradient pulse along z-axis. 1/4 J i2 denotes the evolution 
caused by H s for 1/4 J i2 without pulses. Because of the difference of abun- 
dance of 13 C1 and 13 C2 in the sample and the effects caused by decoupling, 
the equilibrium state of the system is represented as 7 C /* + 7 C / 2 , instead 
of 7c (i^ + I z ), where 7 C is the gyromagnetic ratio of 13 C nucleus. We call 
7(k and 7 C "effective gyromagnetic ratios" for 13 C1 and 13 C2, respectively. 
a = arccos (7c/7c)> where 7^/7^ can be measured by experiment. 

The pulse sequence [f ]*' 2 - ^ - [ir]l< 2 - ^ - [|] 2 transforms | |>i | [> 2 
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into an entangled basis state [13] 



o (0) = I l I 2 - l x l 2 - J 1 / 2 

ys\V) ± x ± x ± z ± z 1 y 1 yl 



(3) 



As soon as the system lies in the entangled state, the decoupling pulses are 
closed. The system and environment evolute under H s+e . By applying a 
hard (nonselective) pulse [vr]*' 2,3 in the middle of a period of evolution time 
t ( shown in Fig.2 ), the effect of chemical shift evolution during this period 
can be cancelled [14] [15]. The system evolutes under H e f represented as 



where the last two terms describe the interaction between the system and 
environment. Only these two terms cause the system to evolute, because 
entangled basis states are the eigenstates of /]J 2 [16]. After t, the environ- 
ment degrees of freedom are introduced into the density matrix of the system 
described as 



p s (t , ID = cos( Vl3 + y? 23 ) (4 1 / 2 - Ill 2 y ) - + 2 S m(y? 13 + ^23) + /;/*)/*, 

(5) 



where ip i3 = 7r J 13 t, and ip 2 3 = ^3?$. By tracing over the environment degrees 
of freedom, the density matrix of the system is obtained, which is represented 



H e f = 27rJ 12 /]/ 2 + 27rJ 2 3/ 2 /f + 27rJ 13 /J/ 3 , 



(4) 
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as 



Ps(t) =3<T \Ps{tjl)\ T>3 +3 <i \Ps{tjl)\ i> 3 = CO S ( V 9 13 +^ 3 )(^4 2 -/ y 1 /, 2 )-/ 2 1 / z 2 , 

(6) 

using Il\ t> 3 = || T>3, and j£| |> 3 = -±| |> 3 . In NMR experiments, Eq.(6) 
is equivalent to the deviation density matrix [13] represented as 



Ps(t) = 



-cos(y?i3 + if 2 ?) 







\ 







(7) 



2, 



y -cos((p 13 + if 23 ) 1 j 

where basis states are arrayed as | |>i I t>2, I T>i I l>2, I |>i I T> 

I |>i I l>2- For convenience, we also use p s (t) in Eq.(7). The function 

cos(ipi3 + ip 23 ) describes the decoherence behavior which can be observed 

in the NMR spectrum through a readout pulse [|]^ which transforms the 



system into the state represented as 
/ 



Psr(t) = 







— I 



icos(ip 13 + if 23 ) cos(ip 13 + if 23 ) 
cos((p 13 + ip 23 ) -icos(ip 13 + ip 23 ) 



-icos(ip 13 + if 23 ) cos((p 13 + ip 23 ) 
cos(ifi 3 + ip 23 ) icos(<p 13 + ip 23 ) 





—2 







(8) 
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Experimental data are taken at controlled temperature (22°C) with a 
Bruker DRX 500 MHz spectrometer. The coupling constants J12 = 103.1Hz, 
J23 = 201.3Hz, and J i3 = 9.23Hz. X H nucleus is again decoupled during 
recording the FID signal. In our experiments, decoupling is the process of 
tracing over the environment degrees of freedom to get the reduced density 
matrix of the system, i.e., Eq.(6) is obtained from Eq.(5). Fig. 3a is the 
carbon spectrum through the readout pulse when t = 3.50ms. The center 
frequencies of 13 C1 and 13 C2 are 124.16 ppm and 117.00 ppm, respectively. 
Fig. 3b shows that the amplitude of the left peak of 13 C1 varies as t. The 
experimental data points can be fitted as function 5.8cos(2nt/T), where T = 
8.72ms. In theory, T = 9.50ms. The error is about 8.2%. It mainly results 
from the imperfection of pulses, and the effect of decoherence which cannot 
be controlled. 

If environment is composed of multiple qubits, the interaction between 
the system and environment is represented as 

iV+2 

Hi = {^J lk l]l k z + 27cJ 2k I 2 z I k z ). (9) 

fc=3 

Under the evolution induced by H iy the system is described as 

iV+2 N+2 

Ps(t) = Illl II C0S (Vlk + V2k) - - I'yly II C0S (Vlk + V2k), (10) 

fc=3 fc=3 
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through tracing over environment degrees of freedom, where ipi k = irJikt, 
and (p2k = nJ2kt- If N — > oo, the coherence usually approaches to very 
fast,if t ^ 0. 

In our experiments, we examine quantum decoherence of a system in 
an entangled state. It is unnecessary to consider the state of environment. 
In fact, we usually cannot describe the state of environment clearly. The 
interaction between the system and environment introduces the environment 
degrees of freedom into the state of the system. Decoherence is the result of 
tracing over such degrees of freedom. 

This work was partly supported by the National Nature Science Foun- 
dation of China. We are also grateful to Professor Shouyong Pei of Beijing 
Normal University for his helpful discussions on the theories in quantum 
mechanics and Miss Jinna Pan for her help during experiments. 
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Figure Captions 

1. The structure of trichloroethylene. The two 13 C nuclei and one 1 H 
nucleus represent qubits 1, 2, and 3, respectively. 

2. The scheme used to cancel the chemical shift evolution, t denotes the 
evolution time. Nonselective 7r pulses are applied in the middle of t. 

3. The experimental results when the two 13 C nuclei in an entangled state 
evolute under if e / (see the text). 1 H nucleus is decoupled during record- 
ing the FID signal. Fig. 3a is the carbon spectrum through a readout 
pulse [f when t = 3.50 ms. The center frequencies of 13 C1 and 13 C2 
are 124.16 ppm and 117.00 ppm, respectively. Fig.3b shows that the 
amplitude of the left peak of 13 CT varies as t. The experimental data 
points can be fitted as function 5.8cos(2nt/T), where T = 8.72ms. 

[Figure 1 about here.] 

[Figure 2 about here.] 

[Figure 3 about here.] 
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